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ABSTRACT 

Laminar  flow  heat  exchangers  which  cool  oil  in  non-inter- 
connected parallel  passages  can  experience  non-uniform  flows  and 
a  reduction  in  the  effective  heat  exchanger  coefficient  in  a  range 
of  Reynolds  number  which  varies  with  tube  length  and  diameter, 
tube  wall  temperature  and  fluid  inlet  temperature.   The  method  of 
predicting  the  reduction  in  effective  heat  transfer  coefficient 
and  the  range  of  Reynolds  number  over  which  these  instabilities 
exist  is  presented  for  a  particular  oil,  Mobil  aviation  oil  120. 
Included,  also,  is  the  prediction  of  the  effect  of  radial  viscosity 
variation  on  the  constant  property  magnitudes  of  friction  and  heat 
transfer  coefficient. 
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INTRODUCTION 

Heat  exchangers  that  operate  in  the  laminar  flow  regime  with  non- 
interconnected  parallel  passages  can  depart  from  predicted  performance  as 
a  result  on  non-uniform  flow  distributions.   Variations  in  the  flow  dis- 
tribution among  the  tubes  can  be  caused  by  the  effect  of  the  fluid  vis- 
cosity/temperature relation,  natural  convection  and  non-uniform  passage 
sizes  (from  design  and  manufacturing  tolerances) .   If  the  exchanger  is 
oriented  vertically  with  the  colder  side  up,  the  effect  of  natural  con- 
vection can  be  eliminated.   Non-uniform  passage  size  can  be  minimized 
through  quality  control.   Therefore,  the  effect  of  the  viscosity/tempera- 
ture relation  on  non-uniform  flow  will  be  analyzed  in  this  paper. 

The  viscosity/temperature  relation  for  the  Mobil  aviation  oil 
120  used  in  this  analysis  is  presented  in  Fig.  1.   Viscosity  is  a  strong 
function  of  temperature,  decreasing  as  temperature  increases.   In  this 

analysis  it  is  assumed  that  other  fluid  properties  such  as  density  and 
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thermal  conductivity  are  uniform  at  p  =  55.0  Ibm/ft   (881.0  kg/m  ), 

k  =  0.08  BTU/hr-ffF   (0.138  W/m'C)  and  cp  =  0.5  BTU/lbm-F   (2.09  x  103 
J/kg*C  ).  This  is  a  reasonable  assumption  for  oils.   Also,  it  is  assumed 
that  pressure  drop  entrance  and  exit  effects  can  be  neglected.   The 
analysis  was  performed  for  flow  in  tubes.   Uniform  wall  temperature  was 
ussed  assuming  the  heat  transfer  coefficient  was  very  much  greater  on 
the  outside  than  on  the  inside  of  the  tube. 

When  cooling  liquids  or  heating  gases,  non-uniform  flow  can 
exist  in  the  flow  passages.   For  the  same  overall  pressure  drop,  two 


discrete  flow  rates  can  exist  in  the  same  heat  exchanger. 

This  problem  has  been  discussed  by  Mueller  [1],  [2]  and  [3] 
and  Rohsenow  [4]. 

EFFECT  OF  RADIAL  VISCOSITY  VARIATION 

The  fully  developed  fluid  and  temperature  distributions  in  a 
tube  are  non-parabolic  due  to  the  viscosity/ temperature  relation.   This 
causes  a  lower  velocity  of  the  colder  fluid  near  the  tube  wall  and  a 
higher  velocity  of  the  hotter  fluid  near  the  centerline  compared  with 
the  constant  property  case  (ie.  uniform  viscosity)  with  the  same  flow 
rate  and  bulk  temperature. 

For  steady  state,  fully  developed  flow  and  viscosity  as  a 
function  of  temperature,  the  momentum  equation  reduces  to: 


dp   13,    du.  ,,  N 

dx   r  dr      dr 


Neglecting  the  effects  of  axial  conduction  and  assuming  constant  heat 
flux  at  the  tube  wall,  the  temperature  distribution  for  fully  developed 
conditions  is : 

rkf7=fe/BcPuTrdr  (2) 

An  expression  is  needed  to  account  for  the  change  in  friction 
factor  and  Nusselt  number  for  the  variable  properties  case  (ie.  variable 
viscosity)  compared  with  the  constant  properties  case  (ie.  uniform 
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viscosity) .   Assume  that  the  effect  can  be  expressed  as  the  ratio  of  the 
wall  to  bulk  viscosities  raised  to  an  exponent  as  follows: 


f     /% 
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%    W 
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Nu 
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Nu 
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..  ■  (k) 


(4) 


f    and  Nu    are  based  on  constant  properties  at  the  bulk 
CP         CP 
temperature. 

To  find  exponents  m  and  n,  Eqs.  (1)  and  (2)  were  solved  simul- 
taneously using  an  iterative  procedure.   Ah  average  fluid  velocity,  wall- 
temperature  and  parabolic  temperature  distribution  were  assumed  initially, 
The  fluid  viscosity  versus  temperature  relation  of  Fig.  1  was  closely 
approximated  by  a  piecewise  curve  fit  using  a  sixth  degree  polynomial. 
Based  on  the  assumed  temperature  distribution,  the  viscosity  y(r)  was 
determined.   Equation  (1)  was  solved  for  the  fluid  velocity  based  on  the 
viscosity  u(r).   The   dp/dx  term  was  adjusted  to  obtain  the  desired 
average  fluid  velocity.   Next,  the  temperature  distribution  from  Eq.  (2) 
was  determined  using  the  calculated  fluid  velocity.   Then,  the  friction 
factor  and  Nusselt  number  were  calculated.   The  viscosity  y(r)  was  de- 
termined for  the  revised  temperature  distribution  and  the  procedure  was 
repeated  until  the  calculated  friction  factors  and  Nusselt  numbers  agreed 
within  less  than  1%. 

The  results  of  this  analysis  showed  that  the  exponents  m  and 


n  varied  as  a  function  of  wall  and  bulk  temperature  as  shown  in 
Figs.  2  and  3.   The  exponents  were  invariant  to  changes  in  the  average 
fluid  velocity. 

SINGLE  TUBE  PRESSURE  DROP 

Having  quantified  the  radial  effects  of  the  viscosity/tempera- 
ture relation  on  friction  factor  and  Nusselt  number  using  the  ratio  of 
the  wall  to  bulk  viscosities  and  the  exponents  m  and  n,  an  analysis 
is  performed  on  a  heat  exchanger  tube.   An  energy  balance  at  the  tube 
wall  for  uniform  wall  temperature  yields  : 


w  c  dT  =  hTT  D(T  -  T  )dx  (5) 

P  o 


or 


TTk   /  hD  \ 
" wc  \ k  / 


o 


Preliminary  calculations  showed  that  the  fluid  flow  was  fully 
developed  a  short  distance  from  the  tube  entrance  and  entrance  effects 
could  be  neglected.   However,  entrance  effects  on  the  temperature  dis- 
tribution had  to  be  taken  into  account  because  of  the  slowly  developing 
temperature  distribution.   The  following  equation  was  suggested  by 
Kays  [5]  to  account  for  variation  in  the  Nusselt  number  in  the  entrance 
region  of  tubes  for  uniform  wall  temperature: 


^_!fp 

Nu=3.66+°-0668^  k 2/3  (7) 

'4  %' 
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Therefore,    the  equation   to   determine   the   temperature   drop   along   the  heat 
exchanger  is   given  by  combining  Eqs .    (4),    (6)    and    (7): 


,      wc 
.„,  IV    \n  /  0.0668—  -r2-  v      , 

'?■)     3.66  +  2x_Jc_   ^_  ix  (8) 


T"To         W/  1+o.04[^^]2/VMC, 

7TX       K 


Equation  (8)  includes  the  radial  effects  and  the  developing  temperature 
distribution  effects.   The  temperature  change  along  the  tube  was  solved 
by  iteration  of  values  of  bulk  viscosity,  local  temperature  and  exponent 
n  .   The  temperature  change  and  bulk  viscosity  along  the  tube  for  various 
flow  rates  is  shown  in  Figs.  4  and  5  for  a  0.25"  diameter  tube  and  tube 
wall  temperature  of  40  F. 

It  can  be  shown  that  the  pressure  drop  for  an  increment  dx 
along  a  tube  is  given  by: 

iE.  =64  1_G2_ 

dx    Re  D  2p  K    ; 

Equation  (9)  is  multiplied  by  Eq.  (3)  to  account  for  radial  effects  on 
the  friction  factor. 


d£     =64_   /Vim    lGi 
dx  R%    \uj       D2p 


In  integral   form  Eq.    (10)   becomes 
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128     f   m   1-m 


The  curves  of  pressure  drop  versus  flow  rate  for  various  lengths 

of  tube  were  calculated  (Fig.  6)  for  D  =  0.25",  T  =  40  F  and  T.  =  300  F. 

It  is  observed  that  AP  goes  through  a  maximum,  then  a  minimum,  followed 

by  a  continuous  rise  as  the  flow  rate  is  increased.  From  Eq.  (11), 

-4 

AP  -  D   ;  therefore  AP  vs.  w  curves  for  other  tube  diameters  would 

look  like  those  of  Fig.  6  but  be  displaced  vertically  proportional 

to  D   . 


MULTI-TUBE  HEAT  EXCHANGER  PERFORMANCE 

These  curves  are  for  single  tubes  but  may  be  used  to  predict 
the  flow  maldistribution  and  resulting  effective  Nusselt  number  for 
multi-tube  heat  exchangers.   This  exchanger  performance  is  determined 
by  how  the  total  flow  is  supplied  -  constant  flow  pump  (gear  pump), 
constant  head  source  or  variable  head-flow  pump  (centrifugal  pump). 

a)   Constant  Flow  Pump 

In  the  case  of  a  heat  exchanger  with  five  tubes  (L  =  4  , 
D  =  0.25"),  non-uniform  flow  can  exist  as  shown  in  Fig.  7.   For  an  in- 
let flow  rate  of  25.  75  lbm/hr  (average  flow  rate  of  5.15  lbm/hr  per 
tube)  and  inlet  temperature  of  300  F  and  wall  temperature  of  40  F, 
four  of  the  tubes  have  a  flow  rate  of  1.26  lbm/hr  with  an  exit  tempe- 
rature of  42.7  F  and  one  tube  has  a  higher  flow  rate  of  20.7  lbm/hr 
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with  an  exit  temperature  of  220.0  F.   The  average  bulk,  temperature  at 
the  exchanger  outlet  is  185.2  F.   Figure  8  shows  a  plot  of  temperature 
along  the  exchanger  for  one  tube  at  higher  flow  and  four  tubes  with 
lower  flow.   The  average  exit  temperature  (Tavg)  is  the  summation  of 
the  individual  tube  exit  temperature  weighted  by  the  percentage  of 
flow  through  each  tube.   This  is  the  temperature  that  would  be  mea- 
sured at  the  heat  exchanger  exit  after  mixing  the  five  streams.   If 
it  were  not  known  that  this  maldistribution  existed  and  it  were  as- 
sumed that  the  flow  was  uniform  in  the  five  tubes  at  one-fifth  of 
the  total  flow,  the  effective  Nusselt  number  would  be  tetermined 
using  Tavg  as  follows : 


L, 

h    .J) 

eff 

u   cc 
eff 

k 

eff  = 

c    W 

P 

n  ttDL 
o 

(12) 


log 


T  -  T. 
o     1 

T  -  T 
i  o    avg/ 


For  this  case  h  _.  =  5.72  BTU/hr-ft2*F  (32.49  W/m2-C)  and  Nu  --.  =  1.49. 
eff  eff 

Also  shown  in  Fig.  8  is  the  temperature  along  the  exchanger  and  the  exit 

temperature  (T   ...   )  if  the  tubes,  in  fact,  had  uniform  flow.   For 
unitorm 

this  case  h '  -  »  10.37  BTU/hr*ft2-F  (58.87  W/m2-C)  and  Nu  --  =  2.70. 

ert  err 

Since  T    >  T   .  _   ,  the  Nusselt  number  (and  heat  exchanger  per- 
avg    uniform 

formance)  is  lower  when  non-uniform  flows  exist. 

An  enlarged  view  of  the  region  of  non-uniform  flows  is  shown 
in  Fig.  9.   On  Fig.  9  curve  (c-d)  at  any  AP  level  represents  the 
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average  flow  rate  with  four  tubes  having  a  flow  rate  on  curve  (u-a)  and 

one  tube  on  curve  (m-1),  the  stable  flow  curves.   Then  5w  ,  =  4w   +  w  „. 

cd     ua    ml 

Curves  (s-f ) ,  (g-h)  and  (o-j)  are  similar  curves  for  flow  splits  of  3-2, 

2-3  and  1-4  tubes  on  curves  (u-a)  and  (m-1) .   As  the  total  flow  rate  is 

increased,  uniform  flow  exists  through  the  peak  at  point (a)  and  on  to 

point  (b).   At  point  (b)  the  pressure  drop  decreases  suddenly  to  point 

(c)  where  one  tube  has  a  higher  flow  rate  while  the  other  four  tubes 

have  lower  flow  rates.   The  total  flow  rate  is:  5w.  =  5w  =  4w  +  w  . 

b     c     u    m 

As  flow  is  further  increased,  flow  increases  in  each  of  the  tubes 
(along  the  pressure  drop  curve)  with  the  average  flow  rate  shown  as 
line  (c-d) .   Point  (d)  is  the  point  depicted  in  Figs.  7  and  8.   As  flow 
is  increased  slightly  above  the  flow  corresponding  to  point  (d) ,  the 
pressure  drop  falls  from  points  (d)  to  (e)  and  now  two  tubes  have 
higher  flow  rates  and  three  tubes  have  lower  flow  rates.   As  the  flow 
is  increased,  the  flow  in  each  tube  increases  with  the  average  flow 
rate  shown  on  the  line  (e-f ) .   This  process  continues  until  uniform 
flow  in  all  five  tubes  occurs  at  point  (k) .   In  the  case  of  infinite 
tubes,  the  average  flow  rate  would  follow  the  line  (a-1)  with  flow  in 
the  tubes  corresponding  to  the  flow  at  either  point  (a)  -  (lower  flow 
rate)  -  (1)  -  (higher  flow  rate).   The  mix  of  tubes  with  higher  and 
lower  flow  rates  is  determined  by  the  average  fluid  flow  rate. 

For  the  case  of  decreasing  flow  rate,  non-uniform  flows  also 
occur.   As  flow  is  decreased  to  point  (m)  the  pressure  drop  decreases 
with  all  tubes  having  uniform  flow.   The  pressure  drop  suddenly  in- 
creases from  points  (m)  to  (n)  at  which  one  tube  has  lower  flow  and 
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four  tubes  have  higher  flow.   As  flow  is  decreased  further,  flow  in  each 
tube  decreases  with  the  average  flow  rate  given  by  line  (n-o) .   Addi- 
tional tubes  "pop-over"  to  lower  flow  at  points  (o),  (q),  (s)  and  (c). 
All  the  tubes  have  the  same  flow  rate  at  point  (b).   In  the  case  of 
infinite  tubes  the  average  flow  rate  would  follow  line  (m-u)  with  flow 
in  the  tubes  corresponding  to  the  flow  at  either  point  (u)  -  lower 
flow  rate  -  or  (m)  -  higher  flow  rate. 

Figure  10  shows  the  effects  of  non-uniform  flow  on  the  Nusselt 
number  as  a  function  of  Reynolds  number  for  the  case  of  five  tubes. 
The  solid  line  is  for  increasing  flow  and  the  dashed  line  is  for  de- 
creasing flow.   The  Nusselt  number  is  reduced  significantly  between 
points  (b)  and  (k)  for  increasing  flow  and  between  points  (m)  and  (b) 
for  decreasing  flow.   These  are  the  regions  of  the  curve  where  non- 
uniform flow  exists. 

For  the  case  of  infinite  tubes,  the  effect  of  non-uniform  flow 
on  exchanger  performance  is  shown  in  Figs.  11  and  12.   As  the  tube 
diameter  is  decreased,  the  region  of  non-uniform  flow  is  extended 
to  higher  Reynolds  numbers  but  the  Nusselt  number  decreases  less. 
Increasing  the  tube  length  has  an  effect  on  exchanger  performance 
similar  to  that  of  decreasing  the  tube  diameter. 

b)   Constant  Head  Source 

For  the  constant  head  source,  as  the  pressure  drop  is  in- 
creased beyond  point  (a)  of  Fig.  9,  the  average  flow  rate  will  in- 
crease to  point  (1)  regardless  of  the  number  of  tubes.   With  con- 
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tinual  increase  in  constant  head,  operation  can  not  be  sustained  at 
flow  rates  between  points  (a)  and  (1) .   For  the  case  of  continual  de- 
creasing constant  head,  all  tubes  remain  operating  along  line   (1)  to 
(in).   Then  the  flow  drops  to  point  (u)  in  all  tubes. 

c)  Variable  Head-Flow  Pump 

For  a  centrifugal  pump  source  with  the  pump  characteristic 
shown  as  lines  (a-b)  and  (c-d)  in  Fig.  13,  non-uniform  flows  can  also 
exist.   The  situation  is  much  the  same  as  for  the  case  of  a  constant 
head  source  except  that  the  average  flow  follows  the  pump  characteristic 
instead  of  the  constant  head  line.  As  pump  speed  is  increased  above 
the  point  (a)  all  tubes  shift  to  operation  at  point  (b) .  The  ex- 
changer can  not  operate  between  points  (a)  and  (b)  at  the  higher 
speed  nor  between  points  (c)  and  (d)  at  the  lower  speed. 

DISCUSSION 

Variations  in  the  tube  wall  and  inlet  temperature  can  affect 

the  values  of  exponents  m  and  n   (Figs.  2  and  3).   Together  with 

the  y  vs.   T  relation  (Fig.  1),  changes  in  exponents  m  and  n 

can  cause  differences  in  the  shape  of  the  AP  vs.   w  curve  determined 

by  Eq.  (11).   The  AP  vs.  w  curves  for  a  single  tube  (L  =  4', 

D  =  0.25")  are  plotted  in  Fig.  14  for  different  combinations  of  wall 

and  inlet  temperature.   The  tube  with  T  =  40°F  and  T.  =  300°F 
r  o  1 

exhibits  a  region  in  which  non-uniform  flow  can  exist.   If  T.   is 
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changed   from  300 °F   to    200 °F,    there  is  only  a  slight   dip   in   the  curve 
and   the  region  of  non-uniform  flow  is   narrower.      If     T       is   changed 

a  O 

from  40°F  to  100°F  with  T.  =  300°F,  no  region  of  non-uniform  flow 
exists. 

Results  were  similar  for  calculations  done  using  another  oil, 
Mobil  automotive  oil  HD-30. 


CONCLUSIONS 

Laminar  flow,  parallel  passage  heat  exchangers  can  experience 
maldistribution  of  flow  and  a  reduction  in  the  effective  heat  transfer 
coefficient  in  a  range  of  Reynolds  number  which  varies  with  length, 
tube  diameter  and  the  levels  of  the  temperature  -  wall  and  inlet. 

For  Mobil  aviation  oil  120  with  the  u  vs.   T  relation  of 
Fig.  1,  the  exponents  m  and  n  in  the  viscosity  correction  ratios, 
Eqs.  (3)  and  (4), are  shown  to  vary  with  both  wall  and  bulk  temperatures, 
Figs.  2  and  3. 
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NOMENCLATURE 

c  specific  heat  at  constant  pressure 
P 

D  tube  diameter 

f  friction  factor  (variable  viscosity) 

fc  friction  factor  (uniform  viscosity  evaluated  at  the  bulk 
P 

temperature) 

G  mass  flow  rate 

h  -,  effective  heat  transfer  coefficient 
ef  f 

k  thermal  conductivity 

L  tube  length 

m  exponent  for  friction  factor  correction 

n  exponent  for  Nusselt  number  correction 

n  number  of  tubes  * 

o 

Nu  ff  effective  Nusselt  number  (variable  viscosity) 

Nu  Nusselt  number  (variable  viscosity) 

Nuc  Nusselt  number  (uniform  viscosity  at  the  bulk  temperature) 

P  pressure  along  the  tube 

AP  tube  pressure  drop 

r  tube  radius,  radial  direction 

Re,  Reynolds  number  evaluated  at  y, 

T  temperature  as  function  of  r  or  x 

T  average  exit  temperature  in  non-uniform  flow  condition 
avg        &         v 

T,  bulk  temperature 

T  wall  temperature 
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T.  inlet  temperature 

T  . ..       average  exit  temperature  for  uniform  flow  condition 
uniform 

u  fluid  velocity  in  x  -  direction 

W  total  flow  rate 

w  individual  tube  flow  rate 

x  direction  along  tube 


Greek  Symbols 

p       fluid  density 

y      viscosity  as  function  of  T 


u       viscosity  evaluated  at  T 
0  o 


y,      viscosity  evaluated  at  T 


b 
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